Temperature dependence of the magnetic Casimir-Polder interaction 



H. Haakh, F. Intravaia, and C. Henkel 

Institut fiir Physik und Astronomie, Universitdt Potsdam, 
Karl-Liebknecht-Strafle 24/25, D-144'^6 Potsdam, Germany 

S. Spagnolo and R. Passante 
CNISM and Dipartimento di Scienze Fisiche ed Astronomiche, 
Universitd degli Studi di Palermo, Via Archirafi 36, 1-90123 Palermo, Italy 



o 

(N 

o 



(N 
> 
cn 
m 

d 

O 



X 



B. Power and F. Sols 

Departamento de Fisica de Materiales, Facultad de Ciencias Fisicas, 
Universidad Complutense de Madrid, Plaza de las Ciencias 1, E-28040 Madrid, 

(Dated: 16 Oct 2009) 

We analyze the magnetic dipole contribution to atom-surface dispersion forces. Unlike its electrical 
counterpart, it involves small transition frequencies that are comparable to thermal energy scales. A 
significant temperature dependence is found near surfaces with a nonzero dc conductivity, leading to 
a strong suppression of the dispersion force at T > 0. We use thermal response theory for the surface 
material and discuss both normal metals and superconductors. The asymptotes of the free energy 
of interaction and of the entropy are calculated analytically over a large range of distances. Near 
a superconductor, the onset of dissipation at the phase transition strongly changes the interaction, 
including a discontinuous entropy. We discuss the similarities with the Casimir interaction between 
two surfaces and suggest that precision measurements of the atom-surface interaction may shed light 
upon open questions around the temperature dependence of dispersion forces between lossy media. 

PACS numbers: 03.70.+k - theory of quantized fields; 34.35.-|-a - interactions of atoms with surfaces; 
42.50.Pq - cavity quantum electrodynamics; 42.50.Nn - quantum optical phenomena in conducting media. 



I. INTRODUCTION 

Ever since the work of Lennard- Jones has the in- 
teraction between atoms and surfaces been of interest in 
many fields of physics, chemistry and technology. The 
seminal work by Casimir and Polder [2| demonstrated 
that the shift in atomic energy levels close to a con- 
ductor is a probe for the quantum fluctuations of the 
electromagnetic field, a key concept of quantum electro- 
dynamics (QED). In this context, a nonzero tempera- 
ture becomes relevant for several aspects of the atom- 
surface interaction: thermally excited motion of the sur- 
face (phonons) and inelastic scattering of atomic beams 
occupation of excited atomic energy levels, and 
enhancement of field fiuctuations due to thermal pho- 
tons 0. The latter aspect is usually associated with 
distances from the surface larger than the thermal wave- 
length Kt = hc/AirkBT, approximately 0.6/im at room 
temperature. The free energy of interaction typically 
shows a change in power law with distance around this 
point: generally, it is enhanced with respect to zero tem- 
perature and becomes proportional to T. This is the 
classical limit where the interaction is mainly entropic in 
character 0-Q. 

Experimental progress in recent years has achieved the 
sensitivity required to detect the small energy shifts that 
occur at distances on the order of At, making use of the 
exquisite control over the motion of atomic beams (beam 



deflection 13, U I , quantum reflection 
of ultracold laser-cooled atoms [13, [l5 |. 



12l |13|) or clouds 
The latter can 



be handled precisely in miniaturized traps implemented 
near solid state surfaces known as atom chips 
These devices use optical or magnetic fields for trapping 
and can hold atomic clouds at distances down to a few mi- 
crons. Here, the atom-surface interaction manifests itself 
typically as a distortion of the trapping potential (loss of 
atoms through tunneling to the surface or change in the 
trap oscillation frequency). Therefore, the design of such 
setups requires exact knowledge of atom-surface inter- 
actions and conversely, theory predictions can be tested 
experimentally with high precision. 

A surprising result of the research on atom chips is the 
importance of magnetic field fiuctuations near the surface 
arising from thermally excited currents in the material 
of the chip (Johnson noise related to Ohmic dissipation) . 
These fluctuations couple to the spin magnetic moment 
of the trapped atoms and are known to provoke the loss 
of atoms from the trap by flipping the sign of the poten- 
tial 1^, 2^ . These losses are a main obstacle for techni- 
cal applications and further down-scaling of atom chips. 
Predictions that superconducting materials reduce the 
spin-flip induced losses significantly have recently been 
backed by trap lifetime measurements 2ll426|. 

In this paper, we address the magnetic dipole contri- 
bution to the atom-surface (Casimir-Polder) interaction 
including nonzero temperature. One would expect this 
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to be a small correction to the electric dipole coupling 
27h32|| because of the smallness of the transition ma- 
trix elements 33|, |3j|. Yet, the strong magnetic mode 
density close to a metallic surface 35-33| and experi- 
mental evidence for magnetic spin flips call for a recon- 
sideration of the magnetic contribution. In addition, the 
thermal occupation of photonic modes is quite relevant 
because magnetic transitions occur at much lower fre- 
quencies than electric ones, leading to a stronger tem- 
perature dependence. Finally, it is well-known that dis- 
persion forces between dielectric and magnetic materials 
can be repulsive [s^ , as has been shown for the magnetic 
Casimir-Polder interaction at T = in Ref. 37|. We were 
thus led to investigate whether at distances beyond At 
the total atom-surface interaction might be reduced due 
to the magnetic contribution. 



In this work, we calculate the magnetic Casimir-Polder 
free energy of interaction at different temperatures and 
consider a few well-known models for the electromag- 
netic response of the surface. Since Ohmic losses are 
crucial for the thermal behavior, it is highly interesting 
to compare both normal metals and superconductors. 
The latter are described here in the frame of the two- 
fluid model and Bardeen-Cooper-Schrieffer (BCS) the- 



ory [3y|, |40[ . We demonstrate that the magnetic atom- 
surface coupling has very peculiar features unknown from 
its electrical equivalent. We find that for normal conduc- 
tors at nonzero temperature, the magnetic dipole contri- 
bution to the interaction is reduced, while it is enhanced 
for superconductors and in certain non-equilibrium sit- 
uations. This resembles the macroscopic Casimir inter- 
action between two dissipative plates, where the correct 
calculation of the force at large distances and nonzero 



temperatures has been the subject of debate [4l| - |43| . 



This article is organized as follows. In Sec. [TTl we 
give a brief review of the formalism used to calculate 
atom-surface interactions. Section IIIII presents the spe- 
cific forms of the response functions and the experimen- 
tal setups they describe. We also give expressions for the 
Green's tensor in different asymptotic regimes of the dis- 
tance between the atom and the surface. The magnetic 
Casimir-Polder free energy and entropy of an atom near 
metallic or superconducting surfaces at zero temperature 
is calculated in Sec. IIVI Section fVl covers the effects at 
nonzero temperature and discusses the dissipative reduc- 
tion in the interaction and questions connected to the 
entropy. Non-thermal (out-of-equilibrium) states of the 
atoms that occur typically in experimental setups are in- 
vestigated in Section IVIIl We conclude summarizing and 
discussing the main results. Further technical details are 
given in the appendices. 



II. ATOM-SURFACE INTERACTION 

Quite a lot of research has been done on the interaction 
between an atom and a surface [1, 27-34 1. It can been 
shown from perturbation theory with respect to the mul- 
tipolar atom-field coupling [27] that the free energy of 
a polarizable particle at nonzero temperature T has the 
following general form (Einstein summation convention) 



T = / duj coth 

2n 



huj 

2fcsT 



lm[l3fju;)n,,{L,uj)]. (1) 



Here, /3 is the (magnetic or electric) polarizability ten- 
sor for the atom in a thermal state of temperature T, and 
is the (corresponding) Green's tensor in the presence 
of the surface, defined in Eq.(|3|) below. In the planar ge- 
ometry we are interested in, the Green's tensor depends 
only on the atom-surface distance L and on frequency. 
It is well known that Eq.(IT]) has the same form for elec- 
tric or magnetic dipole couplings [HI, 131 ; our notation 
is adapted to the magnetic case. A simple and general 
derivation of Eq.([T]) can be given following Refs.plH |4(|. 
The effective interaction potential between a polarizable 
particle and the (magnetic or electric) field B '35[ is given 

by 



{fi-B{ro))T 



(2) 



The expectation value (■■•)t is taken in an equilib- 
rium state of the non-coupled system at temperature T 
and implicitly evaluates symmetrically ordered operator 
products; /x is the (magnetic or electric) dipole operator 
and B the corresponding field operator, evaluated at the 
atom position Tq. The factor 1/2 arises from a coupling 
constant integration (excluding a permanently polarized 
atom). Within first-order perturbation theory, both the 
dipole moment /i(t) and the field B{t) can be split into 
fiuctuating (fl) and induced (in) parts: the fiuctuating 
part describes the intrinsic equilibrium fiuctuation, while 
the induced part arises in perturbation theory from the 
dipole coupling [47|. Eq.JS]) becomes 



(/.-(t).B«(ro,0)T (M«(0-B'"(ro,t))T 



(3) 



Here, we assume the fluctuating parts of the dipole and 
of the fleld to be decorrelated at this order. This assump- 
tion would break down at higher orders of perturbation 
theory. Note that while in Eq.©, the total dipole and 
fleld operators (Heisenberg picture) commute at equal 
times, this is no longer true for their 'in' and 'fl' con- 
stituents in Eq.(j31). The induced quantities are given, in 
frequency space, by the retarded response functions [35| 



/ir(a.)=/3,,(c.)B«(ro,^) 
Br(r,^)=H,,(r,ro,L.)A.«(c.) 



(4) 
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where the frequency dependence aUows for a tempo- 
ral delay. The equilibrium fluctuations follow from the 
fluctuation-dissipation theorem 



where /3° is the state-specific polarizability [2^, [5C 



(^f(..)/i«(c.'))T 



— Siu} — oj') coth I — ^ 
xIm[/3,^(^)] . 



(5) 



(6) 



Combining Eqs-dSHS]), we recover Eq.([T]), setting 
'Hij{ro,ro,uj) = HijiL^uj). One uses the fact that the 
imaginary part of both Green's tensor and polarizability 
tensor are odd in iv (retarded response functions). The 
field correlations are needed at the same position Tq, and 
it is easy to remove the divergent free-space contribution 
(Lamb shift) from T, by keeping in the Green's tensor 
only the reflected part of the field 28| . In a planar geom- 
etry, it follows from symmetry that the result can only 
depend on the dipole-surface distance L. Note that the 
Green's tensor can also depend on temperature through 
the surface reflectivity. In a two-level model for the atom, 
the thermal polarizability 0ij{oj) contains a stronger T- 
dependence because of a Fermi-Dirac-like statistics [i^ , 
see Eq. ([25l) below. 

Eq.([T]) is often expressed in an equivalent form using 
the analyticity of 13^ {i^) and 7i{L,uj) in the upper half 
of the complex frequency plane. Performing a rotation 
onto the imaginary frequency axis yields the so-called 
Matsubara expansion [49| 



(7) 



where ^„ — 2'KnkBT /h are the Matsubara frequencies 
and the prime in the sum indicates that the n — Q term 
must be weighted by a prefactor 1/2. Both /3"^(if) and 
7^(^,1^) are real expressions for ^ > 0. 

If the atom is in a well defined state \a) rather than in 
a thermal mixture, we have the expression of Wylie and 
Sipe [iil 



71=0 
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(9) 



Here, fi^'' ~ (a|/ii |6) are the dipole matrix elements, is 
the frequency of the virtual transition \a) — >■ \b) {uJba < 
for a transition to a state of lower energy). Finally, the 
thermal occupation of photon modes 



n{uj) — (^t 



fiuj/ksT 



(10) 



in the second term in Eq.® is the Bose- Einstein dis- 
tribution. At r = it occurs only for excited states, 
for which n{ujha) — > — 1 for ojba < (see Eqs.(4.3, 4.4) of 
Ref. [2^ ) . The real part of the Green's tensor can be given 
an interpretation from the radiation reaction of a classi- 
cal dipole oscillator 128]. Similarly, this term is practi- 
cally absent for the electric Casimir-Polder interaction 
of ground-state atoms because of the higher transition 
frequencies, « {kg /K) 10^... 10* K. 

In the following, we call the Matsubara sum (first line) 
in Eq.® the non-resonant contribution, and the second 
line the resonant one, because it involves the field re- 
sponse at the atomic transition frequency. 



III. RESPONSE FUNCTIONS 

The formalism presented in the previous section is 
quite general and (3 \TL] could represent either the mag- 
netic or electric polarizability [Green's tensor], respec- 
tively. We now give the specific forms of these quantities 
in the magnetic case, focusing on a planar surface and 
specific trapping scenarios. 



A. Green's tensors and material response 

The Green's tensor for a planar surface can be calcu- 
lated analytically. Let the atom be on the positive z-axis 
at a distance L from a medium occupying the half-space 
below the a;?/-plane. By symmetry, the magnetic Green's 
tensor l-L{L,uj) = 'H(ro,ro,w) is diagonal and invariant 
under rotations in the xy-plane: 



00 

H{L,uj) = ^ / kdk K 
Stt J 



r^^(c.,fc) + 



„TM 



(cj, k) [xx + yy] + 2^r^^(w, k)zz 



-2«L 



(11) 
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where /io is the vacuum permeability and xx^ yy and 
zz are the Cartesian dyadic products. We consider here 
a local isotropic, nonmagnetic bulk medium [/i(a;) — 1], 
so that the Fresnel formulae give the following reflection 
coeflScients in the TE- and TM-polarization (also known 
as s- and p-polarization) [sHj 



r™{u,k) = 



e{uj)K — Kr. 



(12) 



where the propagation constants in vacuum 

and in the medium, respectively (roots with Im k < 
0; Re K > 0) 



\ B ^, 



= JP-e(c.)-, (13) 



and k ~ \k\ is the modulus of the in-plane wavevector. 
Note that the magnetic Green's tensor can be obtained 
from the electric one Q by swapping the reflection coeffi- 
cients pJl] 



n 



9{r 



TE ^ ^TM^ 



(14) 



All information about the optical properties of the sur- 
face is encoded in the dielectric function e{uj). We will use 
four different commonly established descriptions, each of 
which includes Ohmic dissipation in a very characteristic 
way. As it turns out, the magnetic Casimir-Polder in- 
teraction is much more sensitive to dissipation than the 
electric one (see Sec IV Cp . This is due to the fact that 
the resonance frequencies in the magnetic polarizability 
f3{uj) are much lower (see Sec lIII C]) . 

The first model is a Drude metal 13511 



where the order parameter rj{T) follows the Gorter- 
Casimir rule At T = 0, the superconductor coincides 
with the plasma model, as is known from the London the- 
ory of superconductivity [5l|. The plasma model is thus 
the simplest description of a superconductor at zero tem- 
perature rather than a model for a normal metal. More 
involved descriptions of superconductors (including BCS 
theory) also reproduce the plasma behavior at low fre- 
quencies (a; well below the BCS gap) and temperature 
close to absolute zero. The full BCS theory of super- 
conductivity can be applied in this context, too, using 
its optical conductivity ^52r,54i] . as recently discussed in 
Ref.[55|. We shall see below (Sec IV B)) . however, that the 
two-fluid model and BCS theory give very close results 
for realistic choices of the physical parameters. 

Our fourth model takes a look at the peculiar case of 
a very clean metal. Here, rather than by impurity scat- 
tering, dissipation is dominated by electron-electron or 
electron-phonon scattering. In these cases, the dissipa- 
tion rate in the Drude formula ([15]) follows a character- 
istic power law 



7(T) (X T", 



n > 1 



(18) 



at small temperatures and saturates to a constant value 
at high temperatures (Bloch-Griineisen law). It is rea- 
sonable to call this system the perfect crystal model. As 
in a superconductor, dissipation is turned on by temper- 
ature, but in a completely different manner. This can be 
distinguished in the atom-surface interaction potential. 



B. Distance dependence of the Green's tensors 



eDr(w) = 1 



uj{ijj -f 17) 



(15) 



where tOp is the plasma frequency and 7 > is a phe- 
nomenological dissipation rate. This is the simplest 
model for a metal with finite conductivity. If 7 is con- 
stant (independent of temperature), the conductivity can 
be attributed to impurity scattering in the medium. 

The second model is the dissipationless plasma model 
e-p\{ijj): here, one sets 7 = in the right-hand side of 
Ea. ([T5|) . This corresponds to a purely imaginary con- 
ductivity. 

In the context of atom chips, the case of a supercon- 
ductor is particularly interesting because dissipation is 
suppressed as the temperature T drops below the critical 
temperature Tc- We adopt here (third model) a descrip- 
tion in terms of the two-fluid model, a weighted sum of 
a dissipationless supercurrent response (plasma model) 
and a normal conductor response 

e,c(w,r) = v{T)e^M + [1 - 77(T)]eDr(^) , (16) 
4" 



V{T) = 



Q{Tc - T) , 



(17) 



For the Drude model, there are three different regimes 
for the atom-surface distance that are determined by 
physical length scales of the system (see Ref.|32]] for a 
review): the skin depth in the medium. 



zvr V w 



(19) 



where Ap = 2'Kc/ujp is the plasma wavelength, and the 
photon wavelength in vacuum. 



A. 



27rc 



(20) 



Note that e(w) « 2i\1^/ {2-k5^)'^ foi' frequencies w ^ 7 ^ 
LJp (Hagen-Rubens regime). This is the relevant regime 
for the relatively low magnetic resonance frequencies. We 
then have S^j <C Ai^ which leads to the following three 
domains: (i) the sub- skin- depth region, L ^ J^, (ii) the 
non-retarded region, ^ L <C A^j, and (iii) the retarded 
region: <^ L. In zones (i) and (ii), retardation can 
be neglected (van-der-Waals zone), while in zone (iii), it 
leads to a different power law (Casimir-Polder zone) for 
the atom-surface interaction. 
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Since the boundaries of the three distance zones de- 
pend on frequency, the respective length scales differ by 
orders magnitude between the magnetic and the electric 
case. For electric dipole transitions, the Hagen-Rubens 
regime cannot be applied because the resonant photon 
wavelength is much smaller. The role of the skin depth 
is then taken by the plasma wavelength Ap. This implies 
that the Casimir-Polder zone (iii) for the electric dipole 
interaction occurs in a range of distances where magnetic 
retardation is still negligible [zones (i) and (ii)]. 

In the three regimes, different approximations for the 
reflection coefRcients that appear in the Green's func- 
tion pn can be made. We start with the Drude model 
where in the sub-skin-depth zone [T^, we have k 3> 
1/(5 > 1/A and 



„TM 



iuj,k) 



e{uj) + 1 



1 



e(w) 
e(w) + 1 c2fc2 



(21) 



At intermediate distances in the non-retarded zone, the 
wavevector isl/A^fc^l/i5, hence, 



r™(u;,fc) 



.™(c.,fc) 



-1 + i 
l-f-i 



ck 



2 uj 



(22) 



Finally, in the retarded zone we can consider A; <C 1/A <C 
1/(5, so that 



r^^(c.,yfc) « -1 + 



„TM/ 



,k) « 1- 



2 



(23) 



A similar asymptotic analysis can be performed for the 
other model dielectric functions. It turns out that 
Eas. (PTH^5|) can still be used, provided the assumption 
|e(aj)| ^ 1 holds. This is indeed the case for a typical 
atomic magnetic dipole moment and a conducting sur- 
face. 

The asymptotics of the Green's functions that corre- 
spond to these distance regimes are obtained by perform- 
ing the fc- integration in Eg. ljlll) with the above approx- 
imations for the reflection coefhcients. The results are 
collected in Table ID One notes that the zz-component 
is larger by a factor 2 compared to the xx- and yy- 
components. The difference between the normal and 
parallel dipoles can be understood by the method of im- 
ages [sij. Furthermore, the magnetic response for a nor- 
mally conducting metal in the sub-skin-depth regime is 
purely imaginary and scales linearly with the frequency 
w - the reflected magnetic field is generated by induc- 
tion. Only the superconductor or the plasma model can 



reproduce a significant low-frequency magnetic response, 
via the Meifiner-Ochsenfeld effect. In contrast, the elec- 
tric response is strong for all conductors because surface 
charges screen the electric field efhciently. 

The imaginary part of the Green's functions deter- 
mines the local mode density (per frequency) for the mag- 
netic or electric fields [s^. These can be compared di- 
rectly after multiplying by l//io (or Eq), respectively. As 
is discussed in Refs.|36l.l37|. in the sub-skin-depth regime 
near a metallic surface, the field fiuctuations are mainly 
of magnetic nature. This can be traced back to surface 
charge screening that efficiently decouples the interior 
of the metal and the vacuum above. Magnetic fields, 
however, cross the surface much more easily as surface 
currents are absent (except for superconductors). This 
reveals, in the vacuum outside the metal, the thermally 
excited currents from the bulk. 



C. Atomic polarizability 

The magnetic and electric polarizabilities are deter- 
mined by the transition dipole matrix elements and the 
resonance frequencies. We are interested in the retarded 
response function, which for an arbitrary atomic state \a) 
is given by Eq.® above. 

When the atom is in thermal equilibrium, we have to 
sum the polarizability over the states \a) with a Boltz- 
mann weight: 



(24) 



where Z is the partition function. In the limit T — > 0, 
we recover the polarizability for a ground state atom. 
For a two- level system with transition frequency fim, the 
previous expression takes a simple form and can be ex- 
pressed in terms of the ground state polarizability [Eq. 
dH]), where a = g]: 



(3 (cj) = tanh 



2kBT 



13' i^) 



(25) 



Let us now compare the electric and magnetic polariz- 
abilities. The magnetic transition moment among states 
with zero orbital spin scales with fiBds where gs is the 
Lande factor for the electron spin and /is the Bohr mag- 
neton. Electric dipoles are on the order eao with oq the 
Bohr radius. With the estimate that the resonance fre- 
quencies (Om and f2e) determine the relevant range of 
frequencies, we have approximately 



a(o)/£o 1 



(26) 



where ats ~ 1/137 is the fine-structure constant. The 
magnetic interaction is thus expected to be a small cor- 
rection. Conversely, the narrower range of frequencies 
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Subs-kin depth 


Non-retarded 




Retarded 




Drude 


plasma 








Qxx 




MOT 


MO 


MO 


2iLuL 4a;"L" \ 2»^i/, 




16A2L 
1 


327rL3 


327rL3 
1 


^ C (? ) 

f 2ia;L 4a;2L2 \ ai^^/. 




327reoi^ 




327reoi^ 


\ c (9- I 



TABLE I: Magnetic and electric Green's tensors at a planar surface. The other elements have the asymptotes "Hyy = "H-xx, 
Hzz = 2'Hxx, and similarly for Qu. The off-diagonal elements vanish. 



makes it much more sensitive to the influence of temper- 
ature. 

We have seen now that the polarizability of an atom 
takes a positive constant value at low frequency and 
the induced magnetic dipole is parallel to the magnetic 
field (paramagnetism). Let us consider for comparison a 
metallic nanosphere. If its radius R is smaller than the 
penetration dep th and the wavelength, the polarizability 
is given by [35 1 



/3sph(w) 



2tt / Ruj 
15/10 V c 



(27) 



This quantity vanishes at low frequencies and has a neg- 
ative real part (diamagnetism). For a qualitative com- 
parison to an atom one can estimate, e.g., the magnetic 
oscillator strength, defined by the integral over the imag- 
inary part of the polarizability at real frequencies 



j3{uj)dio 



/3sph(w)da; 



n 

27r 



(28) 



^7p) i?^log(^)(29) 



From the Clausius-Mossotti relation follows the electric 
counterpart: 



2^ 

V3 
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(30) 



We find that the nanoparticle has a dominantly elec- 
tric response, similar to an atom, but the ratio of the 
oscillator strengths depends on the material parameters 
and the radius. From the above expressions, we find that 
the absolute value of the nanosphere's magnetic oscilla- 
tor strength is actually smaller than the one of an atom 
if the sphere's radius R ^ Inm. 



D. Optical and magnetic traps 

The resonance frequencies relevant for the magnetic 
Casimir-Polder potential depend on the trapping scheme. 
We focus here on alkali atoms that are typically used 
in ultracold gases and distinguish between optical and 
magnetic traps. 



In an optical trap, we may consider the case that the 
magnetic sublevels are degenerate and subject to the 
same trapping potential (proportional to the intensity of 
a far-detuned laser beam). Magnetic dipole transitions 
can then occur between hyperfine levels whose splitting 
is on the order of r2m/27r « 10®, . . . , 10^° Hz, correspond- 
ing to temperatures of 5,...,500mK (see Appendix [B] 
for more details.) In contrast, electric dipole transitions 
occur in the visible range r2e/27r k, lO^^Hz or ~ 50 000 K. 
If we average over the magnetic sublevels, we get an 
isotropic magnetic polarizability. This allows to write 



Plj = Plo^ij^ so that in Eq.© or © 



Z 



(31) 



The setup we will consider in most of our examples 
is an atom in a magnetic trap. In these traps, one uses 
the interaction of a permanent magnetic dipole with an 
inhomogeneous, static magnetic field B. Let us consider 
for simplicity a spin 1/2 manifold: the Zeeman effect 
then leads to a splitting of the magnetic sublevels by 
the Larmor frequency, VLm = ^Bgs\B\/h in weak fields. 
To give an order of magnitude, 51„i/27r w 280 MHz ss 
{kB/^Trh) 13.5 mK at B = lOmT. Atoms in those mag- 
netic sublevels where AE = —fi ■ B > are weak-field 
seekers, and can be trapped in field minima. The mag- 
netic trap we have in mind is a two-wire trap suspended 
below the surface of an atom chip. Currents in the 
two wires, combined with a static field, create a field 
minimum below the chip surface, with gravity pulling 
the potential minimum into a position where the mag- 
netic field is nonzero and perpendicular to the surface. 
Magnetic dipole transitions are then generated by the 
parallel components fix, l^y of the dipole moment (see 
Appendix IB|) . In this anisotropic scenario, the compo- 
nents of the magnetic polarizability tensor are given by 

PlM) = PlyCO = and PlM) = 0. The rele- 

vant components of the Green's tensor are, therefore. 



Pin 2"^! 



(32) 



We should mention that many experiments do not realize 
a global equilibrium situation, as assumed in Eq.([T]). In 
typical atom chip setups, atoms are laser cooled to /iK 
temperatures or prepared in a well-defined state, while 
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the surface is generally at a much higher temperature, 
even when superconducting. For the description of such 
situations, a more general approach d la Wylie and Sipe 
[Eqs.lIH]) and ([9])] is more suitable, and we discuss the re- 
sults in Sec lVIII Before addressing these, we start with 
thermal equilibrium free energies, however. This may 
be not an unrealistic assumption in spectroscopic ex- 
periments where Casimir-Polder energies are measured 
with atoms near the window of a vapor cell 56|. From 
the theoretical viewpoint, thermal equilibrium provides 
an unambiguous definition of the entropy related to the 
atom-surface interaction. We shall see that this quantity 
shows remarkable features depending on the way dissi- 
pation and conductivity is included in the material re- 
sponse. This closely parallels the issue of the thermal 
correction to the macroscopic Casimir interaction, a sub- 
ject of much interest lately. 

To summarize, in an optical (isotropic) trap, the equi- 
librium Casimir-Polder free energy ([T]) is given by the 
Matsubara sum 

^iso(i,T)= (33) 

OO 

-kB-L Piso(lCn) ^ I 

while in a magnetic (anisotropic) trap, we have 

OO 

J-an(i,T) = -2fcBT^'/3j„(ie„)H,,(L,ie„) . (34) 



IV. ZERO-TEMPERATURE INTERACTION 
POTENTIAL 

The (free) energy vs. distance has been calculated nu- 
merically for an anisotropic magnetic dipole in front of 
a half-space filled with a normal conductor (FiglU top) 
or described by the plasma model (Fig[l] bottom). The 
thick black curves give the zero-temperature result (see 
the caption for parameters). The dashed asymptotes are 
discussed in this section. All energies are normalized to 
the L^^ power law of the non-retarded Casimir-Polder 
energy near a perfectly reflecting surface. The scale fac- 
tor J^pi(l/im), given in the caption, is slightly smaller 
than Eq. ([38| below. These and the following results have 
been obtained from the numerical procedure described in 
Appendix [X] 

The magnetic Casimir-Polder potential is always re- 
pulsive as expected from the interaction between an os- 
cillating magnetic dipole and its image at the conducting 
surface. The sign is also consistent with the macroscopic 
Casimir interaction between a conducting and a perme- 
able surface ('mixed' Dirichlet-von Neumann boundary 
conditions), see e.g. Ref . (ilj . The curves in FigUmake 
manifest the crossovers between the distance regimes in- 
troduced in Sec lIII B] above. The relevant length scales 
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FIG. 1: (Color online) Casimir-Polder free energy vs. distance 
L for an anisotropic magnetic dipole (transition moment par- 
allel to the surface). Drude model (Top) and plasma model 
(Bottom), with plasma frequency ujp/2-K ~ 1.42 x lO^^Hz 
(Ap ~ 210 nm) and 7 = O.Olojp. The transition frequency 
is ^ 480 MHz f» {ks /2Tvh) 0.023 K (A™ = 63 cm). In 

both plots, the free energy scale J-"(l/xm) = 9.79 x 10~^^ J = 
(27r?i)1.48 mHz is the value at T = 0, L = 1/xm obtained in 
the plasma model. Gray dashed lines indicate the asymptotic 
laws at zero temperature Eas. H37H39|) . while colored dashed 
lines include thermal corrections according to Eas. (|44l |45|| . 



are here the skin-depth S„i = 6n,„ , evaluated at the tran- 
sition frequency (for the normal conductor), the plasma 
wavelength Ap = l-Kcjijjp (for the plasma model), and 
the transition wavelength A„j = A(ri„i). The case of 
the superconductor is discussed in Scc lVBI bclow (FigI2|): 
within the two-fluid model adopted here, it is identical to 
the plasma model at zero temperature. The temperature 
dependence interpolates between the Drude and plasma 
case, as discussed in Secs lV Al and IV C[ 

The zero-temperature (black curves) case for a Drude 
model has been stated earlier in Ref. [371]; we give details 
on the asymptotes. Taking the limit T — in Eq.([71) 
recovers the well-known expression 



J"(L,0) = E{V) 



(35) 



In the sub-skin-depth regime L <ti Sm, the distance 
dependence in the anisotropic case p2l) for the Drude 
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model becomes 



•^an(i,0) 



(36) 



where is the magnetic transition dipole matrix el- 

ement, cf. Appendix [BJ This expression is obtained 
by using the sub-skin-depth asymptote of the magnetic 
Green's tensor (first column of Table |I| under the ^- 
integral and cutting the integral off at the border 
of this regime, L ^ (5(^), i.e., ^ « 27c^/(wpL)^. The 
small-distance calculation for the plasma model can be 
done in a similar way. In both the sub-skin-depth and 
non-retarded regimes, the Green's tensor ([TT|) becomes 
independent of ^ (see Table [l| , and the frequency inte- 
gral depends only on the polarizability. Therefore, no 
logarithm appears as in the dissipative case, but 



16A2 



L 



(37) 



In the non-retarded regime (intermediate distances), 
the interaction energy in the Drude model [37] and in 
the plasma model behave alike 



^an(i,0) 



327r 



(38) 



This is calculated as outlined above for the plasma model. 
The energy ([55)1 is identical to the interaction of the mag- 
netic dipole with its image, calculated as for a per- 
fectly conducting surface. Indeed, the power law is 
consistent with the dipole field of a static (image) dipole. 

In the retarded region L ^ A„j, (5„i, Ap (not discussed in 
Ref.js^l), the free energy of the Drude is identical to the 
one of the plasma model. Retardation effects lead to a 
change in the power law with respect to shorter distances, 
identical to the electric Casimir-Polder interaction: 



^an(i,0) 



iMajpMoAm 1 



(39) 



The calculation of this asymptote follows the same lines 
as in the electric dipole case, see Ref . [isj . Comparing 
different transition wavelengths Am (e.g., Zeeman vs hy- 
perflne splitting): the smaller the transition energy, the 
larger the retarded interaction. The numerical data dis- 
played in Fig. [T] agree very well with all three asymptotes. 

Eas. (p8)) and (l39l) illustrate that the magnetic atom- 
surface interaction is reduced relative to the electric one 
by the fine-structure constant a|,, as anticipated earlier 



in Ea. (P5)) . One should bear in mind, of course, that the 
length scales for the cross-overs into the retarded regime 
are very different. A crossing of the non-retarded mag- 
netic and the retarded electric potentials would be ex- 
pected for a distance of order AgO!^ ^ 1 mm, where it 
is clear that both energies are already extremely small. 
In addition, the temperature should be low enough so 



that the thermal wavelength (^i is the first Matsubara 
frequency) 



At = ^- 



fic 



2^1 'i-nksT 



0.18mm K 
T 



(40) 



satisfies A^ ^ 1 mm. Indeed, we shall see in the follow- 
ing section that a nonzero temperature can significantly 
reduce the magnetic Casimir-Polder potential. 



V. CASIMIR-POLDER INTERACTION AT 
NONZERO TEMPERATURE 

In this section, we consider the temperature depen- 
dence of the Casimir-Polder interaction at global equilib- 
rium, in particular using the temperature-dependent po- 
larizability 0^ (PS)) . This provides also a well-defined cal- 
culation of the atom-surface entropy, see Sec lVIl Scenar- 
ios with atoms prepared in specific magnetic sub-levels 
are discussed in Sec lVIIl 

The set of curves in Fig[T] illustrates the strong impact 
of a nonzero temperature for the Drude (normally con- 
ducting) metal: its magnitude is reduced for any distance 
L. In the plasma model (no dissipation), the main effect 
is the emergence of a different long-distance regime: the 
thermal regime L ^ At [Eq. (|40|) ]) where the interac- 
tion becomes stronger than at T = 0. The latter kind 
of behavior could have been expected from the thermal 
occupation of photon modes within the thermal spec- 
trum. The effect in the Drude model is more striking 
and is explained in Sec IV CI below. A significant differ- 
ence with the electric dipole interaction is the fact that 
it is quite common to have temperatures much larger 
than the magnetic resonance energies, fcsT 3> TiU,„i or 
At ^ Am. Thermal effects thus start to play a role al- 
ready in the non-retarded regime, and can be pronounced 
at all distances. 

The usual description of the high-temperature (or Kee- 
som [H^) limit is based on the term n = in the Mat- 
subara sum Q 



J"(L,T-^ oo) 



l3fJ0)H,,iL,0) (41) 



Indeed, the higher terms are proportional to the small 
factor exp(— 2^„L/c) ~ exp(— hL/At) that appears in 
the Green's function T-Lji{L,i£^n). This description is dis- 
cussed in more detail in the following sections. 



A. Plasma model 

In the plasma model and more generally, for all ma- 
terials where the reflection coefficient r'^^{uj,k) goes 
to a nonzero static limit, the magnetic Green's tensor 
HjilLjUj — > 0) is nonzero as well. The leading order po- 
tential in the thermal regime is then given by Eg. (1411) . 
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For the anisotropic polarizability of Eq. ((32)) . and assum- 
ing ksT ^ hflm, the temperature dependence drops out, 
and we find from a glance at Table U 



Mo I Ail 



327rL3 



(42) 



(assuming L 3> Ap). This is identical to the zero- 
temperature result in the non-retarded regime (j38p . as 
can also be seen in FigHJ If the temperature is lower, 
ksT < hQrm but the distance still in the thermal regime, 
the factor tanh(?ir2,„/2fcBr) < Mlm/2kBT in the static 
polarizability reduces the interaction slightly (T = 0.01 K 
in FigUl bottom). 



B. Superconductor 

The atom-superconductor interaction shows a richer 
behavior compared to the plasma model, as illustrated 
in Figl^l At T = 0, it strictly coincides with the plasma 



Si 0.001 - 




FIG. 2: (Color online) Casimir-Polder free energy vs. distance 
for an anisotropic magnetic dipole above a superconductig 
surface. Parameters ujp, 7, fim and normalization F{1 /im) as 
in Figlll critical temperature Tc = 1 K ~ 290 Tm. Tempera- 
tures are T/T^ £ {0, 0.7, 0.9, 0.99, 0.9999, 1.0}. 

model, as it must for the two-fluid description (|17p 
adopted here. The large-distance (thermal) asymptotes 
are the same as in the plasma model for T < Tc- The 
reasoning leading to Eq. (|42l) can be applied here as well: 
the response of the superconducting surface to a static 
magnetic field is characterized by a nonzero value for 
r^^luj 0, k) because of the Meifiner-Ochsenfeld effect. 
Although the superconducting fraction decreases to zero, 
proportional to the product r]{T)ujp, the interaction po- 
tential Eq. p^ stays constant because it does not depend 
on this 'effective plasma frequency'. 

This picture also explains the lowering of the sub-skin- 
depth asymptotes in FiglU from Ea. ([37|) . the Casimir- 
Polder potential is proportional to l/A^ H' r](T){ujp/c)^ . 
This gives scale factors « i, . . . , 2 x 10"'* for the cases 



T = 0.7, . . . , 0.9999 Tc, in quite good agreement with the 
numerical data. 

It is worth mentioning that the full BCS theory can 
give results in close agreement with the simple two-fluid 
model we use here. In FiglSl we show the temperature 
dependence of the Casimir-Polder potential (at fixed dis- 
tance L) for the two cases. We choose here a damping 



parameter 7 = 5 x 10 



in the same order of the zero- 



temperature gap A(0) — 3.5 X 10~'^fiu!p « 1.76fcBTc. 

The BCS calculations have been performed using a re- 
cently developed efficient technique of calculating the op- 
tical conductivity at imaginary frequencies 15511 and using 
an approximative form of the gap equation |59l.[60|. Cal- 
culations over a larger parameter range, but restricted to 
T — 0, have been reported by Skagerstam et al. [34 1. 
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FIG. 3: Casimir-Polder free energy for a superconducting sur- 
face, in the two-fluid model (solid line) and the BCS the- 
ory (dots). Parameters ojp, Q,m as in Fig[l] Scattering rate 
7 = 5 X 10~*ajp, critical temperature Tc = 12K — 500 r,„, 
atom-surface distance L = 1 /im. Energies normalized to 
Jpi(L,T ^ 0) = 1.09 X 10"^" J. 

Going back to FigO note that at T > Tc, the su- 
perconductor jumps to a completely different behavior, 
identical to the Drudc metal. This is expected from the 
two-fluid model (|17p . but also in Mattis-Bardeen theory 
where the gap parameter A(T) vanishes above Tc, and 
the optical conductivity cr(a;, T) coincides with the Drude 
model, see Refs. 52h54|. 



C. Thermal decoupling from a normal conductor 

As mentioned above, the Drude model and the su- 
perconductor around the critical temperature show an 
unusually strong temperature dependence in the mag- 
netic Casimir-Polder potential. The strong suppression 
at large distances (FigslT]and[2|) arises from the fact that 
the Green's tensor Hij{L,i£^) ^ at zero frequency in 
the normal conducting state. The leading order poten- 
tial (|4ip vanishes, and one has to consider the next term 
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^ = ^1 in the Matsubara sum[7l so that the exponentiaUy 
small factor exp(— L/At) governs the thermal (large dis- 
tance) regime. We call this the thermal decoupling of 
the atom from the (normal) metal. This phenomenon 
is related to low-frequency magnetic fields that pene- 
trate the (non-magnetic) surface. Indeed, the vanishing 



0) could have been expected from the 



of Hij{L,u! 

Bohr-van-Leeuwen theorem [6l|, that states that for 
any classical system, the magnetization response to static 
fields must vanish. Both conditions apply here: the ze- 
roth term in the Matsubara series involves static fields, 
and is also known as the classical limit. Indeed, except 
for the material coupling constants, Eq. ipTj) no longer in- 
volves h, while the next Matsubara terms do (via At). 
The Bohr-van-Leeuwen theorem does not apply to a su- 
perconductor whose response is a quantum effect (illus- 
trated, for example, by the macroscopic wave function of 
Ginzburg-Landau theory), and by extension, not to the 
plasma model, as recently discussed by Bimonte [62 1. 

We now calculate the next order in the Matsubara se- 
ries to understand the temperature dependence of the 
Casimir-Polder shift near a metal. For simplicity, we con- 
sider again the limiting case kgT ^ fifl^ which simplifies 
the polarizability to the Keesom form. 



/3an(iCn 



I Ma; 



n> 1. 



(43) 



In the thermal regime, L ^ A^, we use the large-distance 
limit of the Green's tensor (cf. the retarded regime of 
Table HI . The Matsubara frequency a; = i^i then yields 
the mentioned exponential suppression 



a frequency ^„ ^ This leads to an asymptote similar 
to Eg. (1551) . but with the ratio fiS^L/ksT in the argument 
of the logarithm and an additional factor Hflm/kBT. 



VI. ATOM-SURFACE ENTROPY 

It is well-known that at high temperatures where the 
free energy scales linearly in T [Eq. ipT]) ], the dispersion 
interaction is mainly of entropic origin [58|. More pre- 
cisely, the interaction is proportional to the change in 
entropy of the system "atom plus field plus metallic sur- 
face" , as the atom is brought from infinity to a distance L 
from the surface. We calculate in this section the atom- 
surface entropy according to 



SiL,T) = - — 



dT 



(46) 



This entropy definition is unambiguous for the global 
equilibrium setting of the previous Section. 

The behavior seen in the previous figures indicates sig- 
nificantly different entropies for the surface models, with 
a strong dependence on the presence of dissipation (con- 
ductivity) at low frequencies. This parallels the discus- 
sion of the macroscopic Casimir entropy for the disper- 
sion interaction between two plates, a subject of recent 
controversies, where the Drude and plasma models give 
different answers [4i,[63-l63. The results that follow in- 
dicate that the magnetic Casimir-Polder interaction may 
provide an alternative scenario to investigate this point. 



■ expi-L/Ar) 



(44) 



where Am is the magnetic resonance wavelength (cf. 
FigH]). 

At shorter distances, we have to perform the Mat- 
subara summation. In the regime L <C At ^ A™, we 
consider the non-retarded approximation to the Green's 
tensor and make the approximation exp(— nL/Ay) ~ 1- 
The sum over the polarizability /3an(iCn) can then be done 
with the approximation (1431) . and we get 



(45) 



The scaling is in good agreement with FiglT] 

The crossover into the sub-skin-depth regime is now 
temperature-dependent and occurs where the skin depth 
i5(^i) ^ L. This corresponds to a temperature ksTD ~ 
7i7Ap/(27rL)^. The involved frequency is characteristic 
for the diffusive transport of electromagnetic radiation in 
the metal at wavevectors ^ \/L. 

In the sub-skin-depth regime, the leading order ap- 
proximation to the Matsubara sum involves terms up to 



1.0 



0.6 



P 0.4 
0.2 
0.0 



superconductor 
\ perfect crystal 



Dmde i 




1 100 lO-* 



plasma 



1 100 lO* 

T [units of T„t\ 

FIG. 4: Casimir-Polder entropy of an anisotropic magnetic 
dipole with a plasma (solid), Drude-metal (dashed), perfect 
crystal (dotted) and two fluid superconductor (gray). Param- 
eters chosen as in Figs[T] and (2] The unit of entropy AS' is 
defined in Eq. (|47|) . The inset shows a the data for the Drude 
model and its asymptotics obtained from Eqs. (|44I45[) in log- 
arithmic scaling. 

The atom-surface entropy (|46)) is plotted in Fig|4] for 
surfaces of different material. In all models, the entropy 
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vanishes at high temperatures because to leading or- 
der, the free energy (j4ip becomes independent of T, and 
higher orders vanish exponentially with T . (This feature 
is specific to the thermal polarizability of a two-level sys- 
tem.) It is remarkable that the vanishing of the entropy 
happens at temperatures already much smaller than the 
'geometrical scale' UbTl = Tic/AttL [i.e. T/Tl = 2^iL/c 
with ^1 the first Matsubara frequency]. This points to- 
wards another characteristic energy scale in the atom- 
surface system, discussed below. 

One notes in Fig|3]very small values for a plasma and a 
superconductor, two cases where the dc conductivity di- 
verges. The superconductor shows a narrow, pronounced 
entropy peak at T^: we interpret this as the participa- 
tion of the atomic dipole in the phase transition. Indeed, 
the electromagnetic waves near the surface are slightly 
shifted in phase due to the interaction with the mag- 
netic dipole moment. The atom-surface interaction can 
be thought of a sum over all these phase shifts, similar 
to Feynman's interpretation of the Lamb shift. 

In the Drude model, we observe a broad peak at tem- 
peratures where the thermal energy ksT becomes com- 
parable to the photon energy of the magnetic resonance, 
TiVtm = ksTm- Comparable to this scale for our pa- 
rameters is the diffusive energy fcsT/) ^ hjXp/{2TrLy « 
3.3 UV,,n, introduced after Ea. (H5t . We thus attribute the 
atom-surface entropy to the participation of the atom in 
the thermally activated diffusive motion of charges and 
fields below the metal surface ( Johnson- Nyquist noise). 
This motion involves, at the relevant low energies, mainly 
eddy currents whose contribution to the Casimir entropy 
(in the plate-plate geometry) has been recently discussed 
by two of us As T drops below the diffusive scale 

Td, the eddy currents 'freeze' to their ground state and 
the entropy vanishes linearly in T. 

The dotted curve in Fig|3] corresponds to the 'perfect 
crystal' that has not been discussed so far. It gives rise 
to a nonzero atom-surface entropy in the limit T — >■ 
which is an apparent violation of the Nernst heat theo- 
rem (third law of thermodynamics). This has also been 
discussed for the two-plate interaction [HI, [g^, , but 
the entropy defect here has a different sign (it is nega- 
tive for two plates). The sign can be attributed to our 
atomic polarizability being paramagnetic, while metal- 
lic plates show a diamagnetic response. Using the tech- 
nique exposed in Ref.[67[, the limit of the atom-surface 
entropy as T — can be calculated, with the result for 
an anisotropic dipole: 



AS{L) 



A°„(o)hpL(x,^-^o) 



1 



16n hQfn 



(47) 



provides good agreement for T » Tm ■ 

One can argue along the same lines as in Ref.[6^ that 
the Nernst theorem is actually not applicable for this sys- 
tem, since the perfect crystal cannot reach equilibrium 
over any finite time in the limit of vanishing dissipation, 
7 — > 0. The entropy AS" then describes the modifica- 
tion that the atom imposes on the ensemble of field con- 
figurations that are 'frozen' in the perfectly conducting 
material. 

In the two-plate scenario considered in Refs.jiilei] it 
has been shown that not only dissipation but also non- 
locality of the response has strong implications for the 
entropy. In particular, the residual entropy A^* = van- 
ishes, because at very low temperatures, the anomalous 
skin effect and Landau damping take the role of a nonzero 
dissipation rate. Though we have not considered nonlo- 
cality in this work explicitly, one can expect the same 
thing to happen in the magnetic Casimir-Polder interac- 
tion. 

It should be mentioned, that we have also found neg- 
ative values for the atom-surface entropy, albeit very 
small, for temperatures between and and in the 
retarded regime (see Fig 15]) ■ The sign depends on the 
orientation of the dipole and is sensitive to a balance be- 
tween the TE- and TM-polarized parts of the magnetic 
Green's tensor. 
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The second line applies in the non-retarded limit L ^ Ap. 
This expression is used to normalize the data in Fig |3] and 



FIG. 5: Casimir-Polder entropy of a magnetic dipole with 
a surface described in the plasma model, anisotropic and 
isotropic polarizability. Parameters cjp, Q,rn as in Fig[l] atom- 
surface distance L = 1 mm where the thermal effects are bet- 
ter visible. Same entropy scale factor AS as in FigU defined 
in Eq. (|47[) . The anisotropic curve is blown up to become vis- 
ible. 



VII. NON-THERMAL STATES 

We have argued in the last section, that many real- 
istic setups involve non-equilibrium situations. Atom- 
chips are a typical example where two independent phe- 
nomenological temperatures can be introduced for the 
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atom (or a sample of atoms) and the surface. This tem- 
perature gradient is metastable on experimentahy rele- 
vant time scales because of the weak interaction between 
the subsystems. 

We now analyze the case where the surface is described 
by a temperature T, and the atom prepared in a well- 
defined state I a). More complex configurations can be 
studied starting from this simple case. We first consider 
two- level atoms and then multilevel atoms, including hy- 
pcrfinc transitions as they occur for the alkali group. 



A. Two-state atom 

As before, there is a single resonance frequency r2„ 
for the two-level atom. Depending on whether the atom 
is prepared in the ground-state \g) or the excited state 
|e), the sign of the polarizability (0) changes. Referring 
to the first line of Eq.®, /3° does no longer depend on 
temperature. The resonant term of the second line in- 
volves the thermal occupation number that we approx- 
imate by its classical value n{±flrn) ~ '^ksT / {TiVLrn)- 
This is sufficiently accurate at room temperature and 
typical magnetic resonances. In this limit, we obtain a 
simple expression for the magnetic Casimir-Polder free 
energy (anisotropic case, argument 'g' for ground-state 



atom) 

F^L.g.T) « -2kBTY,P'"{^^n)n^AL^£.n) (48) 

ri>l 

+ kBTp^lJ)) {Re [Hxx(i, n^)] - n^L, 0)} . 

The first line is similar to the result in the Drude model 
because of the missing zeroth Matsubara term. From 
this expression, we now discuss the differences between 
the Drude and the plasma model. 

For the Drude model, the resonant contribution involv- 
ing Re'H(i, rim) becomes significant in the non- retarded 
regime. In particular, combined with the non-resonant 
contribution, it changes the sign of the Casimir-Polder 
potential already at short distances, resulting in an at- 
tractive interaction, as soon as T ^ T^, see Figs l6] and 

m 

In contrast, the second line of Eq. p5| nearly van- 
ishes in the plasma model because the magnetic Green's 
function T-LxxiLjUj) is approximately independent of fre- 
quency, at least in the non-retarded regime. We thus get 
a situation where the zeroth Matsubara term is nearly 
removed from the Casimir-Polder potential. The reso- 
nant term still gives the leading order contribution, once 
the expansion of the occupation number is pushed to the 
next order, n(i7m) ~ kBT/Ml„i — \- We then get 

Fll{L,g,T) « -|M,pReH,,(L,r!m) « ^^j^ (49) 



327rL3 

where the last expression applies in the non-retarded 
regime and is identical to the T = case [Eg. (1551) ]. cf. 
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FIG. 6: (Color online) Casimir-Polder free energy vs. distance 
for a two-level atom in its ground state, transition dipole par- 
allel to the surface. The surface is described by a Drude metal 
at different temperatures. Parameters ujp, 7, Q.m as in Fig[T] 
Note the scale factor 1/T to show the classical limit T ^ T. 
The limiting curve at high temperatures (gray dashed line) 
can be inferred from the second line of Ea. l|48p . 
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FIG. 7: Same as Fig[6l but vs. temperature. Distance L = 
1 ^va. The energy scale T{Q) = 2.56 x 10~^*J is the value at 
r = OK. 



FiglSl At larger distances (retarded regime), the differ- 
ence between the Green's functions in the second line of 
Eq. (P5)) is nonzero and becomes the leading term: 



:Ft{L,g,T) 



kuT Mo7r|//^| 
47rL 



hVtr 



47ri 



cos ■ 



An 



47rL 



■ sm ■ 



(50) 



Note that this has a longer range than the power 
law (|3ni): see FigUl This effect is well known from the 
electric-dipole interaction of excited atoms [zOl and con- 
sistent with the classical interpretation (frequency shift 
of an antenna) of the resonant term. 
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FIG. 8: (Color online) Casimir-Polder free energy vs. dis- 
tance for a two-level atom near a plasma surface. Atom in the 
ground state and transition dipole parallel to the surface. Pa- 
rameters Ljp, flm and scale factor J^(l/im) as in Fig[lIbottom) . 
At high temperatures, the curve can be well approximated by 
the classical contribution (thick gray dashed line). The os- 
cillating part of Eq. (|50p not shown here, sets in at distances 

L > Xm. 



B. Trapped rubidium atom 

The atom-surface potential now involves transitions 
|a) \b) to both higher and lower energy levels. Eq.([5]) 
yields the following form of the free energy 



(51) 



n>l,j 
b\2 



+ ^bTY^ j^ {Re [H,,{L,u,,)] - H,,(L,0)} , 



where we assume again that ksT ^ ?i|wfca| which is valid 
in many experiments. The sign of the interaction depends 
on the relative weight of virtual transitions to lower and 
higher energy levels. From Eq.([9]), a virtual level Ei, > Ea 
gives a positive contribution to the polarizability and a 
positive prefactor for the second line in Eq. (j5ip . these 
terms being negative for levels Ef, < Ea. We can inter- 
pret this sign change from the difference between stimu- 
lated emission into the thermal radiation field (for the 
excited atom) and photon absorption (for the ground 
state). Generally speaking, these contributions do not 
cancel each other because the matrix elements lAif P are 
not the same. 

Let us consider the example of *^Rb, prepared in the 
magnetically trappable hyperfine state \a) = \F,mp) = 
|1, —1) of the 5s ground state configuration (Fig|9]). This 
atom has vanishing orbital momentum L = 0, nuclear 
spin / = 3/2, and a single valence electron so that J — 
S = 1/2. The splitting between the hyperfine levels F = 
1 (lower) and F = 2 is rihf/27r « 6.8 GHz, to which the 
Zeeman splitting in the magnetic trap must be added 



with the correct Lande factor. We use the same Larmor 
frequency as before for the two-level atom; because of 
^m/^hi ~ 0.07, we are still in the weak-field regime. 



5s. 



F=2 



F=l 



+2 
+1 
m = 
-1 
-2 



m = 
+1 



FIG. 9: Energy scheme for *^Rb. 

Assuming a quantization axis perpendicular to the sur- 
face, see Sec lIII D] we get an anisotropic polarizability. 
The necessary matrix elements are calculated in Ap- 
pendix|B] Numerical results for the Casimir-Polder inter- 
action according to Eqs.iHl) and ([9]), are shown in Figs.fTOl 
and llll for the Drude and plasma models, respectively. 
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FIG. 10: (Color online) Magnetic Casimir-Polder free energy 
near a Drude metal for a *^Rb atom in a given hyperfine 
state {\F,mF} = |1, —1)). Parameters ujp, 7 as in Fig[l] The 
Larmor frequency (for virtual transitions between neighboring 
Zeeman levels) has the same value ilrn/^-K = (fcs/27r?i) 23 mK 
as before, and the hyperfine splitting is fihf /27r ~ 6.8 GHz ~ 
(A;B/27r;i)0.3K. 



Near a normal conducting surface described by the 
Drude model, the interaction for T = is attractive 
at all distances. We associate the sign reversal (com- 
pared to the absolute ground state considered so far) to 
the coupling to the lower-lying Zeeman levels. At high 
temperatures, the interaction becomes dominated by the 
resonant contribution that grows linearly with T and is 
repulsive. Again, we find the opposite sign as for the 
ground-state atom in FiglHl Thus, the sum of both con- 
tributions leads to a maximum of the free energy at a 
nonzero, T-dependent distance. 

In the plasma model, the potential crosses over globally 
(for all distances) from attractive to repulsive. For all 
practical temperatures, T > Tm, and the interaction will 
be repulsive as shown in FigHT] The results illustrate 
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FIG. 11: (Color online) Same as Fig llOl but for the plasma 
model. The numerical calculations for extremely low temper- 
atures < T ^ Tm, where the potential changes sign, could 
not be performed with sufficient precision. 



that the magnetic dipole interaction of an excited atom 
will be repulsive in all practical realizations. 

To summarize, the strong dependence of the thermal 
correction on dissipation in the surface occurs in both 
non-equilibrium situations considered here: two-level or 
multilevel atoms prepared in a given energy state. The 
magnetic Casimir-Polder potential thus offers an oppor- 
tunity to distinguish between the two models on the basis 
of experimental data taken at low surface temperatures 
and small distances within the possibilities of today's ex- 
periments. 



VIII. CONCLUSION AND DISCUSSION 

We have considered the interaction of a magnetically 
polarizable particle with a metallic or superconducting 
surface including the effects of nonzero temperature and 
out-of-equilibrium situations. Previous work had con- 
sidered mostly the case of electric polarizability, e.g. , 
Ref.fs^l, or was limited to a static magnetic dipole [SSI 
or zero temperature 



34, 37 



The magnetic atom-surface interaction is repulsive 
over a large range of parameters and turns out to be 
highly sensitive to both thermal fluctuations and dis- 
sipation. In this respect, it shows similarities with 
the Casimir interaction between metallic or magneto- 
dielectric plates. 

The T = results of Ref. 37] suggested that the mag- 
netic interaction might be enhanced by raising the tem- 
perature, possibly creating a regime where it dominates 
over the electric contribution. In fact, thermal enhance- 
ment occurs only near a superconductor at distances 
beyond the thermal wavelength, where the material re- 
sponse is governed by the Meissner effect. In normal 
conductors, field penetration prevents such a regime - in 
accordance with the Bohr- Van Leeuwen theorem ~ and 



the Casimir-Polder energy is exponentially suppressed in 
global equilibrium. 

This behavior can be understood qualitatively from 
the competition between attractive and repulsive contri- 
butions to the force. Repulsion arises from the fluctua- 
tions of the magnetic dipole, coupled to its mirror image. 
This is similar to the interaction between electric cur- 
rents (Lenz rule). Field fluctuations, on the other hand, 
produce attractive forces, due to the paramagnetic char- 
acter of the atom polarizability. Both contributions differ 
in their temperature dependence and depend on the state 
of the atom (thermalized, spin polarized). For example, 
attractive forces arise between a ground-state atom and a 
normal conductor, as the temperature scale exceeds the 
magnetic transition energy. Under realistic conditions, 
this flips the sign of the interaction in the regime acces- 
sible to experiments (FiglB]). 

Considering the Casimir-Polder entropy, we found that 
atoms probe the fluctuations in the material: for in- 
stance, at the superconducting phase transition a pro- 
nounced peak appears. In most situations, the entropy 
vanishes at absolute zero in agreement with Nernst's the- 
orem. The only exception we found was the partic- 
ular case of a 'perfect crystal' [Drude dissipation rate 
7(r) oc T", n > 1], already discussed thoroughly in 
the context of the two-plate Casimir interaction 4^, |43| . 
Here, the entropy takes a (positive) nonzero value at zero 
temperature and can be understood as the participation 
of the atom in the disorder entropy of currents frozen be- 
low the surface [g^ . Indeed, the magnetic dipole moment 
mimics the current response of a second plate, except for 
the sign. 

If we compare with an electric dipole, the temperature 
dependence of the Casimir-Polder interaction is relatively 
weak there. This is due to the larger value of the transi- 
tion frequency, so that for realistic temperatures there is 
no difference between the ground-state and the thermal- 
ized polarizability. This was assumed in Ref. . where 
it was also stated that all good conductors behave very 
similarly. In fact, the electric dipole coupling is domi- 
nated by TM-, polarized modes which cannot penetrate 
into the bulk due to screening by surface charges. 

We have shown here that the Casimir-Polder interac- 
tion between a metal and a fluctuating magnetic dipole 
resembles in many respects the Casimir interaction be- 
tween two metallic plates. In both scenarios, the thermal 
dependence is much more pronounced for a normally con- 
ducting surface compared to a description without Ohmic 
losses. The role of dc conductivity in the two-plate sce- 
nario is still an open question at nonzero temperature. 
The atom-surface interaction may thus provide an alter- 
native way to investigate the temperature dependence of 
the Casimir effect, e.g. in atom chip experiments. The 
main challenge is the small value of the interaction en- 
ergy as compared to the electric one. Future measure- 
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ments of the magnetic Casimir-Polder interaction may 
involve high precision spectroscopy on the shift in hy- 
perfine or Zeeman levels. In order to separate effects of 
the magnetic and the electric dipole coupling it will be 
important to find control parameters that affect exclu- 
sively the magnetic contribution. This may include the 
variation in external fields, isotope shifts, and highly po- 
larizable atomic states like Rydberg atoms. 
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Appendix A: Numerics 

The thermal free energies have then been calculated 
numerically in the Matsubara-formalism by summing 
over a thermal spectral free energy density g{i£,) = 
—kBT(3ijHij{i£,,L) at discrete imaginary frequencies. 
Now, the infinite sum is replaced by a finite one plus 
an integral to approximate the remaining infinite partial 
sum 

T{T) = ^ g{zU « E + / Mn)dn . (Al) 

n=0 n=0 •' ^ 

If the upper summation limit N is chosen correctly, the 
remainder is a slowly varying function and the error is 
small, according to the Euler-MacLaurin formula. In all 
systems considered in this work, the exponential in the 
Green's function (jll[) ensures this property. 

Numerical calculations require an automatic and fast 
estimate for the summation limit TV. In a similar scheme 
a fixed number of Matsubara terms [N = 10) has re- 
cently been proposed [6J|. It should be pointed out, that 
since the Matsubara frequencies are linear in T, any such 
scheme with a fixed number of terms breaks down at low 
temperatures, where the integrand g{i£,) has not suffi- 
ciently decayed. A basic criterion for TV is that the par- 
tial sum is sufficiently large, so that the integral is only 



< u < 1 



(A2) 



a small correction 

To avoid the evaluation of the integral in Ea. (|A2[) we 
have used an upper bound 



Ngii^N)/T 



< u 



(A3) 



which exploits the exponential decay of g introduced by 
the Green's tensor. Typically, r ^ L/ At- In our numer- 
ical calculations, we targeted errors u G [10"'^, 10^^] and 
obtained sums over N g [10^, 10^] terms. In some cases 
like BCS theory, the calculation of the optical response 
is not trivial and the remainder integral was neglected 
completely. This requires a small enough value of u and 
yields a systematic numeric error 0{u). 



Appendix B: Magnetic transition matrix elements 

This section gives the magnetic transition matrix ele- 
ments used in the calculations of polarizabilities and free 
energies according to Eqs.dHini) or (PH) . 

The simplest approach used is the spin ^ system (two- 
level atom) with states \ms) = | ± |). The dipole op- 
erator is /X = ^.gTigsS and the quantization axis is z,so 
that 



Sx\ms) 

Sy\ms) 
Sz\ms) 



\\ -rns) , 
inis\ ~ ms) 
ms\ms) ■ 



(Bl) 

(B2) 
(B3) 



According to Eq. (j25p we need only matrix elements 
connecting to the ground-state \g) — | ~ ^), the only 
nonvanishing of which yield 

\{g\S^\e)\' ^\{g\Sy\e)\^ ^\{g\S.\g)\' ^ i (B4) 

Since all polarizabilities ([S]), © or ([M)) are proportional 
to the transition frequency, the Sz transitions connecting 
identical states do not contribute and we obtain the result 
ofEq. ([32ll. 

Let us now calculate the matrix elements for the ^^Rb 
atom prepared in a hyperfine s-state. 



(a|^,|6) = {F,mF\ti,\F',m'p} 



(B5) 



The dipole operator fi = iiBh,{gsS + g^L -f gjl) 
^B^gsS, assuming vanishing orbital momentum and be- 
cause of the smallness of the nuclear Lande g-factor 
gi (X me/rrip. We express the states |a), |6) in the un- 
coupled basis of the spin and nuclear spin momenta with 
the help of the Clcbsch-Gordan coefficients 



\F,mF)= 



(B6) 
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where the action of the components of the spin operator 
on a state is known [Eqs. (|BmB3[) ]. Hence, we find 

•mi ,ras 
nil ,fns 

Mz = 2^ ^sCp,^^,^Cp^^/ . (B9) 
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